Vacuum discharge as a possible source of gamma-ray bursts 
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We propose that spontaneous particle-anti-particle pair creations from 
the discharged vacuum caused by the strong interactions in dense matter 
are major sources of 7-ray bursts. Two neutron star collisions or black 
hole-neutron star mergers at cosmological distance could produce a com- 
pact object with its density exceeding the critical density for pair creations. 
The emitted anti-particles annihilate with corresponding particles at the 
ambient medium. This releases a large amount of energy. We discuss the 
spontaneous pp pair creations within two neutron star collision and estimate 
the exploded energy from pp annihilation processes. The total energy could 
be around 10 51 — 10 53 erg depending on the impact parameter of colliding 
neutron stars. This value fits well into the range of the initial energy of the 
most energetic 7-ray bursts. 



PACS number(s): 98.70.Rz; 52.80.Vp; 26.60,+c 



Gamma-ray bursts (GRBs) were discovered accidentally in the late 1960s by the Vela 
satellites. The discovery was announced in 1973 [TJ. Since then, they have been one of the 
greatest mysteries in high-energy astrophysics for almost 30 years. The situation has im- 
proved dramatically in 1997, when the BeppoSAX satellite discovered X-ray afterglow ||, 
which enabled accurate position determination and the discovery of optical || and radio 
0] afterglows and host galaxies. The distance scale to GRBs was finally unambiguously 
determined: their sources are at cosmological distances 0. In spite of all these recent 
progress, we still do not know what produces GRBs! The nature of the underlying physi- 
cal mechanism that powers these sources remains unclear. The optical identification and 
measurement of redshifts for GRBs allow us to determine their distances and the amount 
of energy that would be radiated in an isotropic explosion. In recent three observations 
(GRB971214 g, 980703 §} and 990123 [§), the total isotropic energy radiated was esti- 
mated to be in excess of 10 53 erg. For GRB990123, the inferred isotropic energy release 
is up to 3.4 x 10 54 erg, or 1.9 M & (where M is the solar mass), which is larger than the 
rest mass of most neutron stars. It has been suggested that the explosion of GRB990123 
is not isotropic, which reduces the energy released in 7-rays alone to be 6 x 10 52 erg 
due to finite beaming angle. However, if one adopts the picture of the fireball internal 
shock model || that random internal collisions among shells produce the highly variable 
7-ray burst emissions, the required initial energy will be raised by a factor of about 100 
since it is argued that only 1% of the energy of the initial explosion can be converted 
into the observed radiation |L(|. Therefore, it appears that the total exploded energy for 
the most energetic bursts is close to or possibly greater than 10 54 erg. It seems to be 
difficult to imagine a source that could provide so much energy. The first and foremost 
open question concerning GRBs is what are the inner engines that power GRBs ? 

On the other hand, the GRB spectrum is nonthermal. In most cases there is a strong 
power law high-energy tail extending to a few GeV. A particular high-energy tail up to 



18 GeV has been reported in GRB940217 |TT| . This nonthermal spectrum provides an 
important clue to the nature of GRBs. 

Various GRB models have been suggested in the literature, see e.g. Refs. [p|, p!2| , p!5| -4l9 
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Among them, the neutron star merger seems to be the most promising candidate. Three- 
dimensional hydrodynamical simulations of the coalescence of binary neutron stars (NS- 
NS) |13|- |15| , ^0|1 , direct collision of two neutron stars |2l|-f23|1 as well as black hole-neutron 
star (BH-NS) merger |24]] have been performed by some authors. The largest energy 
deposition of ~ 10 51 erg by vv annihilation was obtained in the black hole-neutron star 
merger (for NS-NS collision, the total energy is around 10 50 erg |3|). This may account 
for certain low-energy GRBs on the one hand, but it is, on the other hand, still far 
away from the energetic ones mentioned above. However, it should be pointed out that 
in those macroscopic simulations (and almost all GRB fireball models) the effects of 
strong interactions, e.g., the modification of hadron properties in dense matter, many 
body effects, vacuum correlations et al., have been largely neglected except that a nuclear 
equation of state is applied. 

In this Letter we propose an alternative scenario for the source of the most energetic 
7-ray bursts. It is well known that the density is fairly high at the center of neutron 



stars. The central density can be several times nuclear saturation density ||25|| . Further- 
more, superdense matter could be formed at NS-NS/BH-NS mergers and direct NS-NS 
collisions. Three-dimensional hydrodynamical simulations showed that when two neutron 
stars collide with a free-fall velocity, the maximum density of the compressed core can be 
1.4 (off-center collision, the impact parameter b = R, i.e., one neutron star radius) to 1.9 
(head-on collision) times the central density of a single neutron star [^J. At such high 
density, not only the properties of baryons will be modified drastically according to the 
investigation of relativistic mean-field theory (RMF) and relativistic Hartree approach 
(RHA) ||26|| , but also the vacuum, i.e., the lower Dirac sea, might be distorted substan- 
tially |27j since the meson fields, which describe the strong interactions between baryons, 
are very large. At certain densities, when the threshold energy of the "negative-energy 
sea"-nucleons (i.e., the nucleons in the Dirac sea) is larger than the nucleon free mass, 
the nucleon-anti-nucleon pairs can be created spontaneously from the vacuum A 
schematic picture for this phenomena is depicted in Fig. 1. The situation is quite simi- 



lar to the electron-positron pair creations in QED with strong electromagnetic fields 
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The produced anti-nucleons will then annihilate with the nucleons at the ambient medium 
through the NN — > 77 reaction. This yields a large amount of energy and photons. This 
process may happen in addition to the vv annihilation process. The sequential process, 
77 <-* e + e~ , inevitably leads to the creation of a fireball. The dynamical expansion of the 
fireball will radiate the observed 7-rays through the nonthermal processes in shocks ||. 
In the following, we shall estimate whether enough energy is available within this scenario 
to satisfy the requirement of a source of energetic GRBs. 

We start from the Lagrangian density for nucleons interacting through the exchange 
of mesons 

+ \mluj^ - Ir^R/ 4 " + ^m 2 p -R, ■ R^ 



where the usual notation is used as given in the literature [26 . Based on this Lagrangian, 
we have developed a relativistic Hartree approach including vacuum contributions which 
describe the properties of nucleons and anti-nucleons in nuclear matter and finite nu- 
clei quite successfully [^9|]. The parameters of the model are fitted to the ground state 
properties of spherical nuclei. The RHAO set of parameters gives ^{M-n / m^) 2 = 229.67, 
g^iMj^ / m^) 2 = 146.31, g 2 (M N /m p ) 2 = 151.90. It leads to the nuclear matter saturation 
density p = 0.1513 /m~ 3 (0.1484 - 0.1854 /m" 3 ) with a binding energy E bind = -17.39 
MeV (-16 ±1 MeV) and a bulk symmetry energy a sym = 40.4 MeV (33.2 MeV). The cor- 
responding empirical values are given in parentheses. The model can be further applied 
to the neutron-proton-electron (n-p-e) system under the beta equilibrium and the charge 
neutrality conditions which is in particular important for the neutron star. The positive 
energy of the nucleons in the Fermi sea E + and the negative energy of the nucleons in the 
Dirac sea E_ can be written as 

1/2 1 



E_ = - 



I [k 2 + (M N - g a af] V2 + g^ + -gpT 0j Ro,o} , (2) 
I \k 2 + (Mat - g CT cr) 2 ] V2 - g^o + ^g P T i?o,o j • (3) 



Here a, uo and i?o,o are the mean values of the scalar field, the time-like component of the 
vector field, and the time-like isospin 3-component of the vector-isovector field in neutron 
star matter, respectively. They are obtained by solving the non-linear equations of the 
meson fields including vacuum contributions under the constraints of charge neutrality 
and general equilibrium. The energy of anti-nucleons E + is just the negative of EL, 
i.e., E + = —E- [p9fl . By setting k = in Eqs. (2) and (3), one gets the energies of 
nucleons and anti-nucleons at zero momentum. The critical density pc for nucleon-anti- 
nucleon pair creation is reached when E_ = Mjy. The results are given in Fig. 2 where the 
single-particle energies of the positive-energy nucleon and the negative-energy nucleon are 
plotted as a function of density. Due to the effects of the p- meson field, pc = 6.1 p for 
pp pair creation and 7.5 po for nn pair creation. At the same time, we have calculated the 
equation of state (EOS) of neutron star matter. The structures and properties of neutron 
stars can be obtained by applying the equation of state to solve the Oppenheimer-Volkoff 
equation |JT|. The maximum mass of stars turns out to be M max = 2.44 M & , and the 
corresponding radius R = 12.75 km and the central density p cen = 5.0 p . The p cen is 
smaller than the critical density pc- That means that the spontaneous NN pair creation 
does not happen for a single neutron star within the model employed. 

We consider the following case of neutron star collision: Two identical neutron stars 
with p cen = 4.5 po (with the current EOS, it is related to M = 2.43 M Q and R = 13.0 km) 
collide with each other with a free-fall velocity. The impact parameter b stays between 
and R, which determines the factor of density enhancement. We assume that a compact 
object of average density 7.2 p is created in the reaction zone. The radius of the compact 
object is assumed to be r = 1 km (case A) or r = 3 km (case B) depending on the values 
of b. Since for a single neutron star with p cen = 4.5 po the density at r = 1 km is 4.46 po 
and at r = 3 km is 4.18 po, in case A the density is enhanced during neutron star collision 
by a factor around 1.6 while in case B around 1.7. In both cases the pp pair creation 
will happen while the contributions of the nn pair creation is negligible (it contributes 
at higher density but does not affect our discussions). We define a Dirac momentum hp 
which describes the negative-energy nucleons occupying the eigenstates of the Dirac sea 
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from the uppermost level (the lowest-energy antiparticle level) to the negative continuum 
(see, Fig. 1), i.e., E- = —M^ in Eq. (3). At the critical density for pp pair creation 
Pc? = 6.1 po, the Dirac momentum fc£ = 11.28 /m _1 ; and at p = 7.2 p , ho = 12.45 
fm' 1 . We further define a momentum p max at 2£_ = Mjy, which turns out to be 

Pmax = (gu^O - ^g p r i?o,0 + ga^ - 2M N ^j (g^Uo - ^g p T R 0fi - g a o^j . (4) 

Based on the semi-classical phase-space assumption we then estimate the number of the 
pp pairs whose energies are larger than the nucleon free mass at p = 7.2 p as 



N pair = -vrr 3 x ^ = 2.147r 3 x 10 54 . (5) 
1 3 3vr 2 



3 . , Pmax _ o 1 A 7^,3 w i n54 



Before expansion, this compact object remains high density. Let us check whether most 
of the pp pairs can be created spontaneously. The rates for the iViV pair production 
per unit surface area and unit time, dN pa i r /dSdt, has been calculated in Ref. for 
compressed matter. In the case of p = 7 po and tunnel distance d = 1 fm, the rate turns 
out to be 2.68 x 10~ 2 fm~ 3 . For case B with r = 3 km, the time needed to emit the 
available pp pairs is t — 1.9 x 10 19 fm = 6.3 x 10~ 5 s, which is smaller than the typical 
dynamical scale of NS-NS collision r ~ 10~ 3 s. Thus, we have enough time to produce 
proton-anti-proton pairs spontaneously. The produced protons stay in the atmosphere 
due to gravitational force. However, at that time holes (anti-protons) are still in bound 
states due to potentials they feel (a small fraction may be transported into the negative 
continuum). The above process may happen before the compact object expands (we are 
discussing a microscopic procedure in a macroscopic phenomenon). Then the compact 
object expands and the potentials in the Dirac sea fall down. Those anti-particles (holes) 
in bound states are pushed into the lower continuum and thus escape. They annihilate 
with the protons in the atmosphere or in the surrounding objects and release a large 
sum of energy. If one assumes that 80% of the produced anti-protons annihilate with 
protons in the surrounding medium and the released energy is 2 GeV at each event (at 
the moment it's not very clear how many anti-protons in the Dirac-sea can escape through 
the lower continuum. This is a problem, which should be investigated more closely), the 



total exploded energy E to t turns out to be 5.5 x 10 51 erg and 1.5 x 10 53 erg for cases A 
and B, respectively. As mentioned before, the efficiency to transfer the initial energy to 



the observed radiation is only 1% |T0[. It seems to be necessary to adopt the picture of 
beaming explosion for the most energetic 7-ray bursts. 

Some discussion is now appropriate. Neutron star collisions have repeatedly been 
suggested in the literature as possible sources of 7-ray bursts H32lj33f , powered either by 



annihilation or by highly relativistic shocks. In Ref. Ruffert and Janka claimed 
that a 7-ray burst powered by neutrino emission from colliding neutron stars is ruled 
out. Here we propose a new scenario caused by the strong interactions in dense matter. 
A large number of anti-particles may be created from the vacuum when the density 
is higher than the critical density for spontaneous particle-anti-particle pair creation. 
Such high density can be reached during the NS-NS collisions, BH-NS mergers, or even 
NS-NS mergers when the merged binary neutron stars have large maximum densities. 
Some of the produced anti-particles can be ejected from the reaction zone due to violent 
dynamics. They may be the novel source of low-energy cosmic-ray anti-particles which is 
currently an exciting topic in modern astrophysics p3 |. Most of them will annihilate with 
the corresponding particles at the ambient medium, and thus release a large amount of 
energy. As a first step we have discussed the pp pair creation in two neutron star collision 
scenarios because its critical density is lower than that of other baryons. Our calculations 
show that the exploded energy satisfies the requirement for the initial energy of the 
energetic GRBs observed up to now. The variation of the released energies of different 
GRBs can be attributed to the different impact parameters of colliding neutron stars. 
Those anti-protons, although produced spontaneously, annihilate during the dynamical 
procedure with random probability in collisions with protons. Furthermore, the anti- 
protons annihilating later might be accelerated by the photons produced by the nearby 
pp pair annihilations taking place earlier. This leads to the high-energy anti-protons and, 
consequently, the high-energy photons. Some of them may escape from the fireball before 
being distorted by the medium. Those escaping high-energy photons may constitute the 
observed high-energy tail of 7-ray bursts. This has to be pursued further theoretically. 
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In summary, we have proposed a new scenario of vacuum discharge due to strong 
interactions in dense matter as a possible source of 7-ray bursts. Based on the meson 
field theoretical model we have estimated the exploded energy E tot ~ 10 51 — 10 53 erg within 
two neutron star collisions, which fits well into the range of the initial energy necessary 
for most energetic 7-ray bursts. For a more quantitative study, one needs to introduce 



hyperon degrees of freedom [£5[] and even quark degree of freedom f35||36j if one assumes 
that the center of neutron star is in quark phase. Here we have mainly discussed NS-NS 
collisions. In fact, the proposed scenario may happen more frequently for BH-NS mergers 
since the production rate for BH-NS binaries is ~ 10~ 4 per yr per galaxy |37| which is 
much larger than the rate of direct NS-NS collisions (for an estimation of collision rate 
in dense cluster of neutron stars, see Ref. In this case one might obtain a even 



higher explosion energy reaching the value of 10 54 erg. A relativistic dynamical model 
like relativistic fluid dynamics incorporating meson fields is highly desirable to simulate 
NS-NS collisions, NS-NS/BH-NS mergers. At the end, we would like to mention that a 
similar process may happen in nucleus-nucleus collisions as discussed in the introduction of 
Ref. |29[ where a dynamical production of anti-matter clusters due to the variation of the 



time-dependent meson fields has been suggested. We propose to study the photon and 
anti-proton spectra in ultra-relativistic heavy-ion collisions which may provide us with 
information of structure of discharged vacuum. Works on these aspects are presently 
underway. 
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FIG. 2. The single-particle energies of the positive-energy nucleon and the negative-energy 
nucleon in neutron star matter. 
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